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Introduction: Various models of the reflectance 
from particulate surfaces [e.g. 1-2] are used for inter- 
pretation of remote sensing data of solar system objects 
[e.g. 3-4], These models rely upon the real (n) and 
imaginary (X) refractive indices of the materials. Such 
values are limited for commonly encountered silicates 
at visual and near-infrared wavelengths (X, 0.4-5 pm). 
Availability of optical constants for candidate materials 
allows more thorough modeling of the observations 
obtained by Earth-based telescopes and space craft. 

Two approaches for determining the absorption co- 
efficient (a=27tX/X) from reflectance measurements of 
particulates have been described; one relies upon 
Kubelka-Munk theory [5] and the other Hapke theory 
[6-7]. Both have been applied to estimate a and k for 
various materials [e.g. 6-10]. Neither enables determi- 
nation of the wavelength dependence of n, n=f(k). 
Thus, a mechanism providing this ability is desirable. 

Using Hapke-theory to estimate k from reflectance 
measurements requires two additional quantities be 
known or assumed: 1) n=f(k) and 2) d. the sample par- 
ticle diameter. Typically n is assumed constant (c) or 
modestly varying with X [9-10]; referred to here as n 0 . 
Assuming n 0 , at each X an estimate of k is used to cal- 
culate the reflectance and is iteratively adjusted until 
the difference between the model and measured reflec- 
tance is minimized. The estimated k's (X/) are the final 
results, and this concludes the typical analysis [9-10]. 

Analytical Approach: The approach described 
here combines Hapke theory with subtractive Kramers- 
Kronig (sK-K) analysis [11]. Hapke theory is used to 
determine X, and the next step applies a sK-K analysis 
to yield n h where ni=f(k). These steps are repeated 
until no significant change in n and k occurs. 

Experimental: The montmorillonite used here 

(SWy-1) is from the Clay Minerals Society and is 
originally from Wyoming. The <38- and 38-45 pm 
grain size fractions were separated via dry sieving. 
Reflectance spectra were measured with incident and 
reflected beams at 8° and 6°, respectively. 

Results: Using d of each sieve fraction, derived via 
SEM (3.38- and 3.03-pm), and assuming n=1.405 [12], 
X/ s were derived (figure la). Although the k values 
agree well, there is a difference of -100 between these 
and the results of [12]. 

Reevaluation of the SEM image analyses suggest 
c/’s are biased to smaller values as large grains were 
excluded because they either had small grains adhering 
to them, or were aggregates of small grains. 


An approach in derivation of k that initially relies 
upon the median particle size of a sieve fraction has 
been described [10]. Addressing the observed discrep- 
ancy at the longer X’s, the analyses were repeated using 
the median values (19- and 41.5 pm). The resulting k 
values are shown in Figure la. Although these are a 
factor -2-3 different from each other, they are closer to 
the data of [12], 

[10] also described adjusting median particle c/’s to 
force X’s derived from different sieve fractions into 
agreement. Using this approach c/’s of 26 pm for both 
sieve fractions produces k values in agreement with 
those of [12, figure lb]. In the iterative approach de- 
scribed above, d of 26 pm is used. 

Figure lc shows the difference between successive 
values derived for n (An). Going from n=c to n=/(X), 
i.e. n r n 0 , the greatest An is produced, -0.015 (gray 
dashed line). Subsequent iterations, n 2 -«; and n 3 -n 2 , 
show a decreasing An until it is zero. Similar trends 
are seen in Figure Id for the difference in derived k 
values (AX). The procedure quickly reaches a solution 
and further iterations are unnecessary. 

One additional effort was motivated by the desire to 
estimate the optical constants of montmorillonite for 
the broadest wavelength range. It involved combining 
the X’s derived here with those of [12] for the same 
sample (figure If). These X values were used in the sK- 
K analysis to provide n=f(X) (figure le). 

Comparison with previous data: [8] reported a 
for one Wyoming montmorillonite (sample #25). To 
compare a’s figure 3 of [8] was digitized independ- 
ently six times. A linear interpolation was used to cre- 
ate: 1) a common X scale, from 0.4 to 2.5 pm with a 
sampling of 0.01 pm; and 2) associated a’s. From 
these interpolated values, the average and standard 
deviation were calculated (figure 2). 

It is only feasible to provide a semi-quantitative 
comparison of the a’s derived by the two different 
techniques. This is due to the difference in the spectral 
resolution of the instruments used to obtain the data of 
[8] (see [5]) and the measurements of SWy-1. 

Inspection of Figure 2 shows agreement is good for 
the highest a’s within the region of X overlap of the 
two data sets, with values of about 7- and 15-cm' 1 near 
1.4- and 1.9-pm, respectively. In X regions outside 
these higher values, a’s agree to a factor of -2, except 
at the longest wavelengths, where the measurements of 
[8] are only qualitative. The local maximum at -2.2 pm 
in the SWy-1 data corresponds to a hydroxylated sili- 
cate spectral feature [13]. 
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Figure 1. (a) k values of SWy-1 derived using: 1) d 
determined via SEM (solid lines); 2) median of each 
sieve fraction (long dashed lines), for the large (L) and 
small (S) sieve fractions, and values from [12] (short 
dashed line), (b) k values derived with d= 26 |im, for 
large (L) and small (S) sieve fractions, are comparable 
to the values of [12]. (c) An for successive iterations. L 
and R refer to the ordinate axis, (d) same as c, but for 
A k. (e) n via sK-K analysis of k shown in panel f. (f) 
k’s derived and combined with those of [ 12]. 



Wavelength, iim 

Figure 2. Comparisons of the absorption coeffi- 
cients derived in this study for montmorillonite (gray 
points and line) to those determined by [8] (black 
points with error bars). The dashed arrow indicates 
estimated values shown by [8]. 

The discrepancy between the wavelength positions of 
the local a maxima at ~1.9 |im may be due to the dif- 
ference in the intrinsic instrumental resolutions. Alter- 
natively, it may be due to differences in the composi- 
tion or hydration state of the samples. Cationic substi- 
tution and hydration state influence the position, width, 
and asymmetry of the montmorillonite 1.9 lim band 
[14-15], 

Summary: Procedures described here provide the 
ability to estimate the wavelength dependence of both 
n and k from reflectance spectra in the visible and near- 
infrared. The resulting a’s compare favorably with 
those previously reported for montmorillonite [8]. 
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